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Abstract 

$y\ We will in this note show that it is possible to diagonalise the Lund Fragmenta- 

,-G ■ tion Model. We show that the basic original result, the Lund Area law, can be 

factorised into a product of transition operators, each describing the production 
of a single particle and the two adjacent breakup points (vertex positions) of the 
string field. The transition operator has a discrete spectrum of (orthonormal) 

C^ ■ eigenfunctions, describing the vertex positions (which in a dual way corresponds 

to the momentum transfers between the produced particles) and discrete eigen- 
values, which only depend upon the particle produced. The eigenfunctions turn 
out to be the well-known two-dimensional harmonic oscillator functions and the 
eigenvalues are the analytic continuations of these functions to time-like val- 
ues (corresponding to the particle mass). In this way all observables in the 
model can be expressed in terms of analytical formulas. In this note only the 
1 + 1-dimensional version of the model is treated but we end with remarks on 
the extensions to gluonic radiation, transverse momentum generation etc, to be 
performed in future papers. 
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1 Introduction 

The Lund Fragmentation Model is built upon a few very general assumptions: 
there is a string-like force field between the coloured constituents, there is causal- 
ity and Lorentz covariance, the production of the particles can be described in 
terms of a stochastical process and the process will obey a saturation hypothe- 
sis. Using semi-classical probability considerations we are then led, |jj.g], to a 
unique stochastical process for the breakup of the force field into the final state 
hadrons. In Section g, we will provide a set of necessary formulas to describe 
the dynamical developments along the surface, spanned by the string field. The 
major result is that the probability to reach a particular (exclusive) final state 
is given by the phase space of the state multiplied by a negative exponential 
of the area spanned before the string decays ("the Lund area-law"). In gen- 
eral the model has been used as it is implemented in the well-known Monte 
Carlo simulation program Jetset, H. This means, on the one hand, that it 
is possible to take into account a large amount of kinematical complications, in 
particular from the decay of the primary produced resonances. On the other 
hand, in order to make the simulation programs time-effective, it is necessary to 
introduce routines that make the process rather difficult to follow. In particular, 
it is difficult to disentangle the major dynamical features of the model from the 
many necessary numerical compromises in the simulation program. 

In this note, we will show how to diagonalise the basic stochastical process, 
i.e. how to define a complete set of eigenfunctions and eigenvalues describing 
the process on the local level. In this way we can provide analytical formulas 
for all possible correlations between the observables in the process. 

In Section 0, we will factorise the Lund area law in a somewhat different 
way than the ordinary. It is done by defining an operator that describes the 
transition from one production point to the next in the process. It can be written 
as an integral operator describing the probability to go from one breakup point 
along the string field to the next thereby producing a particular hadron. We 
show how to diagonalise the operator in terms of its eigenfunctions (which in 
a very neat way corresponds to two-dimensional harmonic oscillator functions) 
and calculate its eigenvalues, that are closely related to the Lund fragmentation 
function. 

These eigenvalues turn out to be the major building stones in all the model 
correlations and in Section |J, we show a set of properties of the eigenvalues. It 
turns out that the eigenvalues form a discrete set, and that they correspond to 
analytical continuations of the harmonic oscillator eigenfunctions from space- 
like to time-like regions. We show how to use the factorisation properties of the 
model to provide a set of useful relations for the products of the eigenvalues. 
We also exhibit the relationship to a field theory in a two-dimensional Euclidian 
space that will be further pursued in future publications. 

We will in this note be satisfied to treat only the simplest case corresponding 
to the 1 + 1-dimensional dynamics of the Lund Model. We will, however, end 
with an outlook on future work, in particular on the effects of gluon emission 
and transverse momentum properties of the hadronisation process. 



2 The Lund Fragmentation Model 




Figure 1: A high-energy string breakup of an original qoqo -pair, created at a 
single space-time point. 



The massless relativistic string is in the Lund Model used as a model for the 
colour force fields with colour-3 quarks (q) and colour-3 antiquarks (q) at the 
endpoints. Finally, the colour-8 gluons (g) are in the Lund Model interpreted 
as internal excitations on the string. We will, in this note, treat all the (qq)- 
particles as massless and moving along the lightcones. However, the final result 
is independent of this assumption, cf. Q2|, (massive q and q would in a semi- 
classical scenario move along hyperbolas with the lightcones as asymptotes). 
Further, we will as examples of the formalism consider only e + e~-annihilation 
reactions and refrain other processes as well as gluonic bremsstrahlung to future 
work. 

Then an original q g -pair is assumed to be created at a single space-time 
point and start to go apart thereby stretching the string field in between them 
(cf. Fig pi) . The field will break up into new pairs at the vertices Xj = (xj+ ,Xj-) 
(we use lightcone coordinates) and a q = qj-\ will together with a q = q~j from 
the adjacent vertex form a final state hadron with the energy- momentum pj. 

In this way we have introduced a convenient ordering in the form of rank: the 
first rank particle is formed by (q qi), the second rank by ((71(72) etc. It is also 
possible to introduce a rank-ordering from the q side, i.e. along the opposite 
lightcone starting with (q q n -i) in an n-particle final state. The dynamical 
results should of course be independent of the ordering. Actually, it is easy to 
convince oneself that it is necessary that all the vertices must be space-like with 
respect to each other. One finds that the energy-momentum of the j-th particle 



is given by pj — k(xj-i + — Xj+, Xj 



-i_) (here k ~ 1 GeV/fm is the string 



constant and we will for simplicity put it equal to unity) . As the vector pj must 
be time- like (with squared mass equal to m3 = Pj) we conclude that the two 
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Figure 2: The production process described in terms of momentum transfers in 
a chain. 



adjacent vertices are space-like. 

There is immediately a second conclusion. If we define the vectors qj = 
(xj+, — Xj_) then we obtain that pj = qj—i — qj. This means that the Lund 
fragmentation process also can be described by means of a ladder graph as in 
Fig. H Thus the energy-momentum transfers between the particles along the 
ladder is in a dual relationship to the production vertices in the description in 
Fig. |l[ We will use this relationship in Section g to derive the results of this 
paper. 

The derivation of the Lund Model formulas is based upon these observations 
and a final assumption: even if the energy of the original pair (in the cms 
denoted conventionally W = y/s) will increase without limit the distribution of 



the proper times of the production vertices Tj 



Xj+Xj- will stay finite. 



In terms of the momentum transfers this "saturation assumption" evidently 
corresponds to (one of) the ordinary assumptions behind Gribov's Reggeon 
theory, that the momentum transfers stay finite in this limit. 

To see the details, we will concentrate on two adjacent vertices in the center 
of the process with the coordinates Xj and Xj-\ such that the above-mentioned 



hadron with pj is produced in between, cf. Fig. p 
the coordinates z + = 1 — Xj + /xj_i + and Z- = 1 



It is convenient to introduce 
Xj-i—/xj- that are Lorentz 
invariants and will have the range < z± < 1 independent of the other variables. 
We also describe the vertices by the hyperbolic coordinates (Tg, yi ), £ = j — 1, j 
and note that due to Lorentz covariance the process can only depend upon the 
T's (i.e. the proper times squared) and the relative hyperbolic angles 5yj = 
yj-i — yj (note that the <5y's will be fixed by the mass requirements). 

We will then consider the breakup vertex Xj—i to be the last in a long 
row of production points along the positive lightcone and using the saturation 




Figure 3: Two adjacent vertices with coordinates Xj and Xj—i (in the center of 
the process) and a hadron with energy-momentum pj produced in between. 



assumption we expect that the probability distribution is 

H(xj—i)dxj-i+dXj—i- = H(Y j_i)dT j-idyj-i 



(1) 



In the second line we have made use of Lorentz invariance to claim that the func- 
tion H only depends upon r\,-_i. The production of the particle pj is then given 
by another "step" along the positive lightcone with the probability f(z + )dz + to 
take the fraction z + of the remainder. Then the combined probability is given 

by 



H(Tj-i)dTj_idyj-if(z + )dz + 



(2) 



In the same way we may consider the production of the particle as the last in 
a long row of steps along the negative lightcone, firstly arriving at the vertex 
Xj with the probability H and then taking another step along the negative 
lightcone. In this way we obtain the joint probability 

ff(r,)dr j dy J /(z_)dz_ (3) 

The basic assumption in the Lund Model is then that the two probability dis- 
tributions in Eqs. (^ and (J3|) are equal and this defines in a unique way the 
distributions H and /,@,[| 



Hj(T) = CjT^ exp(-&r) 

fi-ijiz) = Wj-i^l - zpz?'- 1 -"'- 1 exp(-6m 2 /z) 



(4) 



The parameters cij (with the notation for fj-i,j meaning that the hadron with 
mass m is produced in a step from the point j — 1 to the point j) may be 
different for different vertices (e.g. spin- and/or flavour dependent) but the 
parameter b should be the same, i.e. it must correspond to a general colour 
dynamical property. (Speculations on its origin can be found in e.g. B)- In the 




Figure 4: The two areas used in Eq. (ph. The large area is the one spanned 
below the first meeting point of the two constituents from the adjacent vertices. 

phenomenological applications of the Lund Model there has (besides the first 
particle in a heavy quark jet according to a suggestion by Bowler, Qj) been no 
use for more than a common a-value. We will in general in this paper treat this 
simpler case and only when it is useful exhibit the differences to the general 
case. Finally, the parameters Cj and Nj—ij are normalisation constants. 
The joint probability distribution H(T)f(z) can then be written as 

H(T)f{z) ex CN[(1 - z)Y\ a exp(-6(r + m 2 / ' z) = (area) a exp(-b(Area)) (5) 

where the two areas, the large and the small one, are shown in Fig. ^. Evidently, 
the opposite production direction will produce the same result and the areas play 
therefore due their simple factorisation properties (just as they do in general for 
gauge field theories) a fundamental role in the process. (If there are different 
values of the a-parameter the result is similar with different areas represented, 
one typical of each vertex) . 

We will next consider the probability to produce a rank-connected iV-particle 
set {pj}, for definiteness along the positive lightcone starting at the turning 
point of the original q (cf. Fig. ||). The first rank particle will then take 
a fraction z\ of the total lightcone energy momentum, the second will take a 
fraction Z2 of what is left, i.e. (1 — z\) etc. The observable fractions are then 
Ci = Zi, C2 = 2 2 (1 - CO, C3 = z 3 (l - Ci - C2) = zs(l - z i)(! ^ z i) etc. In this 
way we obtain in easily understood notations, H,0]: 

N J? 
l[N j ^(l-z J reM-bm 2 j /z j ) = 

1 z ? 

N dt N 

II^T 1 ( 1 -EO) a exp(-^ + r)) = 
1 « 1 

dz ._ ,„ , , ,_ . , , -n- ,, duj 



da — (1 - z) a cxp(-6s(l - z)/z) TT Ni —?- x 

z • y u . 

N N 

cxp(-bA)S(l - Y^ Uj)6(s - Y m 2 j/uj) = dP ext dP int (6) 




Figure 5: An N-particle cluster with notation as explained in the text. 

In the second line, we have introduced the variables Q defined above and in 
the third the common variable z — J^i (j an d finally rescaled the fractions 
into Uj — Cj/z. We have also introduced the total area (according to Fig. ||) 
A to t — A + T and the total cms energy ^/s of the TV-particle cluster. We note 
that the "final" vertex proper time squared is T = (1 — z)s/ z. 

The final result is then that the probability distribution can be factorised into 
two parts. One of them (dP ex t) (note that it is independent of the multiplicity 
N) is the probability to make a cluster of mass -y/i and lightcone fraction z. 
We note the close similarity to the fragmentation function for a single particle 
/ in Eq.(|J). The other one (dPi nt ) is the probability that the cluster will decay 
into just these particular particles with the fractional energy momenta Uj in the 
cluster cms. 

The distribution dP ext for the general case when several values of a occurs 
will only depend upon the first and the last a-values: 



dP ext =ds—z a °( — 

z z 



-) a -exp(-bT) 



(7) 



The distribution dP ext can be used to study the convergence of the "sat- 
uration assumption" on the distribution H . One finds, H that for squared 
cms-energies s larger than a few times the inverse of the parameter b there is 
an exponential convergence. 

The distribution dPi nt can be reformulated using that duj/uj is equivalent 



to d pS(p 



and that the two delta functions in the same way can be 



written in terms of the particle energy momenta (pj 
Ptot — \Ptot+ > s j Ptot+) 



(u j P to t+,m 2 /u : jPtot+)), 



6(1 



1 



n n 

1 1 



Ptot) 



(8) 



Putting this together we obtain the Lund Model area law: 

n n 

dP int = J] NjtPpjStfj - m 2 3 )S{Y, Pj ~ Ptot) cxp(-6A) (9) 

1 1 

In the case of different a-parameter values we obtain an extra factor for each j: 
u- with (Sa)j = dj—i — a,j, i.e. only the differences of the adjacent a-values 
occur. (It is useful to note that u+jU-j = m 2 / s in order to see that the formula 
is symmetric between the forward and backward lightcones). 

3 The Transition Operator and its Eigenfunc- 
tions 

We will in this section rearrange the Lund Model area law , cf. Eq. (g), in 
another form i.e. as a product of a set of step operators taking us from one 
vertex to the next thereby producing a particle in between. We will after that 
show that this transition operator has a well-defined set of eigenfunctions with 
discrete eigenvalues. 




Figure 6: The area A subdivided into triangular regions as described in the text. 

To do that, we note that the area A can according to Fig. @ be subdivided 
into (hyperbolic) triangular regions (each with an extra "tip" corresponding to 
half of the squared mass; we neglect them for the moment and note that they 
can be included in the particle production constants Nj). The size of these 
regions are given by (SA)j = WTj~iTj\8inh(5yj)\ (we are using the notations 
from the earlier section) and one finds by direct calculation that 



{5A) 3 = yJxiTj^T^-mj)^ 

X(a, 6, c) = a 2 + b 2 + c 2 - 2ab - 2bc - 2ac (10) 



If we introduce the lightcone fraction z according to Fig. o, then we can write 

r,-_i = (l-z)(T j + m 2 j /z) 

(SA), = (zT, + m 2 /z)/2 (11) 

We also note the vectors q.j = (xj+, — Xj_) (with —q 2 = Tj) that fulfil pj — 
qj—i — qj, thereby "solving" the energy- momentum conservation conditions in 
Eq.([|). If we introduce them instead of the particle vectors {pj}, then we find 
the Jacobian 



d 2 P ,S( P 2 - m 2 ) = dTj/y/xtTj-uTj,-™}) (12) 

We can consequently subdivide the whole process according to the arealaw as 

n 

dPint = Y\_ K(T j ^i,T j ,m 2 )dT : j 
1 

exp(-6/2JA(r~UV^|)) 

KiT^T^m 2 ) = Nj V (13) 

A useful representation of the kernel function K is (it is easily obtained from 
the considerations above) 

K(T^ 1 ,r j ,m 2 ) = 
r 1 dzexpi-bizTj + m 2 Jz)/2) 

/ p^_^ iULA 5{T ._ x - {1 - z){T . + m y z)) ( i4) 

JQ z 

(To be precise the result in Eq.(O) must be supplemented by boundary condi- 
tions but we will neglect them because in this paper we will only be interested 
in results outside the fragmentation regions). It is useful to consider the eigen- 
functions of the transition operator K. For simplicity we will introduce the 
dimensionless variables Tg — ► bTg = xg and m 2 —^ bra 2 = yi and consider the 
solution to the equations 

Kg n {x) = K(x,x',y)g n (x')dx' (15) 



The surprising and very gratifying result we obtain is that the functions g n are 
well-known in mathematical analysis. We are going to call them the Laguerre 
functions (noting that the Laguerre polynomials L„ are orthonormal in the 
measure dxexp(-x)) 

g n (x) = L„(x)exp(-x/2) (16) 

They are orthonormal on the positive real axis < x < 00 in the measure dx 
(we use the notations from |H) 



oo 



dxg n (x)g m (x) = 5„ 
'o 

J29n(x)g n (x')=S(x-x') (17) 

n 



Further, the Laguerre functions are the eigenfunctions of the two-dimensional 
harmonic oscillator corresponding to angular momentum equal to zero. In fact, 
it is easy to prove that the eigenfunctions g n will fulfil the following equation 
because of the well-known properties of the Laguerre polynomials L n 

(- A +bQ 2 )g n (bQ 2 ) = 2(2n + l)g n (bQ 2 ) (18) 

In the second line, we have considered the two-component vector Q with the 
scalar product Q 2 = Q\ ± Qq. The differential operator is correspondingly 
defined as A = ^ ± 9 \ e ^ the equation is valid both for Euclidian metric 
and for space-like directions in two-dimensional Minkowski space. 

We note in particular, that the eigenfunctions are independent of the mass 
m of the particle produced in between, this mass-dependence comes solely in 
the eigenvalues 

A„ (y) = N exp(y/2) f dz/z exp(-y/z)(l - z) n (19) 

JO 

Here we have kept to the definition of the hyperbolic triangle, without the 
"tip" , in the kernel function K (cf. the remarks before the Eq.([lC|)). We will be 
concerned with the properties of the eigenvalues in the next section but we note 
at this point their close relationship to the Lund fragmentation function / (for 
equal a- values) in Eq.([|). It is also obvious from Eq.([l^) that the eigenvalues 
will be discrete and decrease quickly with n. 

To obtain these results from Eq. ([15]) , we make use of the representation of 
the kernel function K in Eq.rtlJ) and find the following necessary (and sufficient) 
requirements on the Laguerre polynomials: 

X n L n (x) = dz/zexp(-y/z)L n ((l-z)(x + y/z)) (20) 

Jo 

It is easy to see that for a polynomial of the n:th degree the eigenvalues A„ will 
have to fulfil Eq. (|19|) . However, to prove the general result in Eq.(p0|) we have 
expanded both sides in powers, performed for the m:th term inside the integral 
< £ < m partial integrations (to get rid of the powers of y) and then gathered 
the powers in x. We feel that there must be a simpler way but we have not 
found it yet. 

Given these results, we note that a theorem attributed to Mercer (private 
information, |6|) provides the following representation for the transition operator 
K: 

OO 

K{x,x',y) = ^2g n (x)\ n (y)g n (x') (21) 

71=0 

In order to check the convergence properties of Eq.(Gfi), we show in Fig. M the 
results for the ratio of the left-hand to the right-hand side of the equation. It 
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Figure 7: The ratio R= e * p( M*** £J»/ ^ =0 g n (x)X n (y)g n (x') for different 

X*(x,x ,—y) 

values of p,x and x' when y=0.5. 



is evidently in general only necessary to keep a few terms to obtain a good 
approximation. 

Due to the orthonormality of the Laguerre functions, it is further immedi- 
ately obvious that while Eq. ( pl|) represents the distribution after a single parti- 
cle production between x' and x the result for the production of N particles in 
between them is given by 



K N (x,x') = ^2 9n(x)(X n (y)) N g n (x') 



(22) 



n=0 



In the next section, we will show how to provide a formula for a fixed invariant 
mass-square s and/or a hxed lightcone fraction z (thereby completely defining 
the relationships) between the points labelled by x' and x. 



4 The Properties of the Eigenvalues 

We will in the following discussion for simplicity put the normalisation constant 
N equal to unity but we will insert it in the end formulas. Then the eigenvalues 
X n {y) (with y = bm 2 ) defined in Eq.(|l9|) will have the following property: 



exp(y/2)A n (y) 



1 (fa (i_, ) »55(-( 1 -*W*) 



f dzJ2(i-z) m+n L m (y) = Yl 

J m m 



z 

L m (y) 
m + n + 1 



(23) 
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We have here in going from the second to the third line in Eq. ( |23| ) made use 
of the generating function for the Laguerre polynomials (0): 

Vr u,\ 9 « - v- (-y^ k sr z "' +fc ( m + fc ) ! _ c M-yz/(± - z)] . 

2^L n {y)z _^-^-^ ^ - ^7) (24) 

n k m v ' 

In the second line, we have introduced the series expansion for the Laguerre 
polynomials, rearranged it by changing the original index n — > m + k and then 
summed up firstly a negative binomial and then an exponential series. 

Consequently, the eigenvalue for n can be written as a series in the eigen- 
functions g m and we may from this representation immediately conclude (using 
the differential equation in Eq.(fL8h) that the eigenvalues will fulfil 



(- A +bp 2 )X n (bp 2 ) = -2(2n + l)X n (bp 2 ) + 46{bp 2 ) (25) 

Here we have introduced the timelike vector p with length equal to the mass 
(and the differential operator A is defined in terms of p:s components). We 
have also used the result in Eq.(|17]) noting that g m (0) = 1 for all values of m. 
We conclude that also the eigenvalues are governed by the harmonic oscillator 
equation and that they correspond to a particular analytic continuation of these 
functions from spacelike to timelike vectors. 

Actually, the eigenvalues X n (y) are solutions to the (degenerate) hypergeo- 
metric differential equation and in conventional notations, B in terms of Whit- 
taker functions we have X n (y) — n\W- n -x/2fl{y)/^/y. Using either the dif- 
ferential equation in Eq.(j25|) or the formulas in B we obtain another useful 
representation that bears out these analyticity properties: 

^) = f*fc|Mi« (26 ) 

Jo t + y 

Next we will consider the correlation coefficients for the case when we produce 
N particles in between the vertices denoted x and x' in Eq. (E2n. Using the 
same procedure as in connection with the derivation of the iV-particle cluster 
in Eq.(|6j), we can immediately write 

A^(y)= I dsR N (s)X n (bs) 

r N N 

Rn(s) = / \{N d 2 p 3 5{p] - m 2 )S(Y,P 3 ~ Ptot)eM-bA) 
■* i i 

X n (bs)=exp(bs/2)X n (bs) (27) 

The quantity X n (bs) is the probability (in the nth harmonic oscillator state) to 
produce a cluster with the energy s. It is the integral over all z values of dP ext 
in Eq. (ph (with the Lund parameter a exchanged for n). In the same way 
Rn(s) is related to the integral over dP int in Eq. (ph, i.e. it is the phase space 
integral (including the area law) of the N particles. We have brought back the 
normalisation constants Nj in the expression for Rn- 



11 



N particles with a fixed squared mass s in between: 

^g n {x)\ n {bs)R N {s)g n {x') (28) 

n 

We will next turn to the properties of the phase space integrals Rn but before 
that we make the following observation. The model has very simple factorisation 
properties both in the energy-momentum fractions and, as we have shown in this 
note, in the energy-momentum transfers (or in the dual language, in the vertex 
positions). Consequently, it is in the same way as in Eq.(|2q) possible to pick 
out any other variables, like particular energy- momentum fractions somewhere 
"in between" and reformulate the remaining correlation coefficients accordingly. 
(It is of course necessary to define the scaling variables properly). As every 
possible observable is either related to the energy momentum transfers or to 
the energy-momentum of the observed particles we have in this way a complete 
analytical description of the process. 

We will show two particularly simple and useful properties of the phase space 
integrals Rn . They fulfil a set of iterative integral equations and there are very 
simple formulas for the analytic function TZn which is defined by 

n N{u) ^f ds Mfl (29 ) 

J s + u 

To see these properties we note that if we "pick out" the first rank particle 
from the iV-particle cluster that defines Rn we obtain, cf. Eqs. (g) and ( p7[ ) 
Nidui/uiexp(—bm 2 /ui). The remaining (N — 1) particles will give the same 
contribution but the energy is reduced to si = {Ptot^Pi) 2 = (l — ui)(s — m 2 /ui). 
We have consequently the integral equation 

Rn(s)= 1 Ml i?jv-i((l - ui)(s - m 2 / Ul )) exp(-6m 2 /m) (30) 

Jo Ul 

We note the similarity to the original integral equations for the harmonic os- 
cillator functions g n (x), but there are two major differences. The first is the 
change of sign in front of the term m 2 ju\ in the integral (reflecting the fact that 
we are going from the space-like vectors qj to the timelike vector Ptot)- The 
second difference is that the argument in the functions Rn and Rn~i do not 
have the same range. It is evident that the threshold for producing TV particles 
is SN.thresh — N 2 m 2 which is larger than the threshold for N — 1 particles. 

There is, however, another relation which will make it possible to calculate 
the functions Rn analytically at least as a perturbation series. In order to see 
that, we consider the following sum and use the results from Eq.(p7[): 

Y,(K(y)) N L N (bu) = 

n 

/ dsR N {s) / — exp-bs{l-z)/z^2(l-z) n L n (bu) = 
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f dz 
cIsRn(s) / — exp— (bs + bu)(\ — z)/z = Tt^{u)/b (31) 

We have once again used the generating function in Eq. ( p4| ) and performed the 
z integral. In this way we have a representation of the analytic function 1Z along 
the positive real u-axis. It is , however, necessary to extend this function to the 
negative real axis in order to obtain the properties of Rn- We note in passing 
that the first line in Eq.(pl|) is (besides a factor exp(— bu/2) the contribution for 
the case when we would start out at the lightcone x = and consider N steps 
to the point x 1 = bu (g„(0) = 1 for values of n). This result can be used in a 
perturbation theory by noting that for large values of n we have the following 
behaviour of our functions: 

g n (x) =* J (2v^) and \ n {y) ~ ^2/^K (2^) (32) 

We can then in this approximation write 

Ef d 2 v . 
9n ~" / ~2~ exp ( z2u W 

n J n 

^^^7(2^)3/2 ¥Ty (33) 

We have then defined two-dimensional Euclidian vectors with ft 2 = u and with 
v 2 — » n. Therefore the whole expression can approximately be written 



C xp(-bu/2)K(u) = J2(K(y)) N 9n(bi 



d 2 tj f d 2 v cxp(i2v(fl — J2 tj)) 



N 

n^^^+y) 



2 (£ + v). 2ir 



N , 2 

TT °L^ S (jl - Y &)/4 (34) 

In this way we have exhibited the analytical function exp(—bu/2)lZN(u) for 
values of u > as the contribution from a simple expression obtainable in a two- 
dimensional Euclidian field theory. (Note that we use the Laguerre functions g n 
and not the Laguerre polynomials L n in our approximations). 

We note that the approximation corresponds to the large n-limit of the 
harmonic oscillator function, i.e. we are far away from the ground state and 
consequently "the motion" behaves as almost free oscillations. 

It is possible inside the same formalism to take the neglected terms in the 
approximation into account as further contributions in the model. We can use 
the present results to show that the convergence radius in an expansion of 1Z n 
around u = is given by N 2 m 2 (just as expected). But further terms in the 
expansion are necessary in order to obtain the precise threshold behaviour. Wc 
will present such results in future publications. 
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5 Concluding Remarks 

Due to its simple factorisation properties, the Lund Fragmentation Model can as 
we have shown in here be diagonalised in terms of harmonic oscillator functions. 
We have, up to now, only treated the 1 + 1-dimcnsional version of the model 
but we will, in future publications, continue the work into the 3 + 1-dimensional 
real world. 

Transverse momentum is in the fragmentation process of the Lund Model 
produced via a tunneling mechanism, leading to a gaussian spectrum. In the 
simplest version of the model, there are no correlations between the transverse 
momentum generated at one vertex and at the next but the experimental data 
show, M, that such correlations occur at least in the production of the light 
pions. A mechanism, with strong similarities to the Ornstein-Uhlenbeck process 
for the velocity distribution of a Brownian motion particle, has been proposed 
and succesfully applied to the data, g . 

As this process again is of a gaussian character, we may use very similar 
methods to diagonalise it, cf. also [9J. There is, however, a small subtlety. If 
the transverse momentum is firstly generated and afterwards the string-field 
used to provide longitudinal momentum (as in Jetset) then the transverse 
mass is used instead of the ordinary mass. This would provide a particular 
correlation between transverse and longitudinal motion. We will come back in 
a later publication to a general investigation of the dynamics. 

Gluon radiation is in the Lund Model treated in terms of internal excitations 
of the string field and this will lead to a bent string surface, H. The fragmen- 
tation of states containing one or more gluons has been introduced into the 
Lund Model by Sjostrand, JL0J, using a particular generalisation of the process 
described above. This process is implemented into Jetset and has been very 
succesful to describe the experimental data. The method we have introduced 
above for a flat string surface can be almost directly applied (there is a mi- 
nor change, that we feel may have some implications for the description of the 
fastest particles in a gluon jet) to the method introduced by Sjostrand. 

Bosc-Einstcin correlations has been introduced into the Lund Model by in- 
terpreting the Lund Area Law as stemming from the square of a quantum me- 
chanical transition matrix element, |llf| , [[12|. This means that the Area Law for 
production of two or more identical particles, obtains a weight factor depend- 
ing upon area differences. One problem in this respect is that with n identical 
particles the weight factor will obtain contributions from the n! different per- 
mutations of the particles. The mathematicians call these problems exponential 
in n and they are very time-consuming (although we were in |12] able to bring 
them down considerably). As the area differences arc directly expressible in 
terms of the variables we have discussed above we have some hope to be able 
to reexpress the full result by means of our formalism. 

Finally, we have obtained a new set of tools to study the energy dependence 
of quantities like Rn{s) and even more interestingly the sum over all the mul- 
tiplicities R(s) = J^n R- n - This dependence will necessarily be of the kind s a . 
This is known from before, 0] , H , but the power a, that corresponds to a "Regge 
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intercept" in quark scattering, in accordance with Gribov's Reggeon theory, will 
in this way be accessible for analytical treatment. Evidently the corresponding 
power obtained for a "gluon fragmentation process" will have some meaning for 
the soft Pomeron intercept. 
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